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Abstract 

We develop a QCD sum rule analysis of the form factor F^*^*-,^o(^q^ ^Q"^) in the region 
where virtuality of one of the spacelike photons is small <^ 1 GeV^ while another is 
large: Q^> 1 GeV^. We construct the operator product expansion suitable for this kinematic 
situation and obtain a QCD sum rule for F^*^*^o (0, Q^). Our results confirm expectation 
that the momentum transfer dependence of (0, Q^) is close to interpolation between 

its = value fixed by the axial anomaly and Q^^ pQCD behaviour for large Q^. Our 
approach, in contrast to pQCD, does not require additional assumptions about the shape of 
the pion distribution amplitude (/37r(x). The absolute value of the l/Q^ term obtained in this 
paper favours (/^^(x) close to the asymptotic form ip'^^{x) = 6/^x(l — x). 
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1. Introduction. 



The transition 7*7* 7r° of two virtual photons 7* into a neutral pion provides an exceptional 
opportunity to test QCD predictions for exclusive processes. In the lowest order of perturbative 
QCD, its asymptotic behaviour is due to the subprocess Y ill) + 1* (12) q{xp) + q{xp) with x (x) 
being the fraction of the pion momentum p carried by the quark produced at the qi (52) photon 
vertex (see Fig. la). The relevant diagram is similar to the handbag diagram for deep inelastic 
scattering, with the main difference that one should use the pion distribution amplitude (Ptt{x) 
instead of parton densities. For large Q^, the perturbative QCD prediction is given by 



' ' 3 Jo xQ^ + xq^ 3 jo xQ^ 3Q^ 



{Q^ = — g|, q^ = —q\ and our convention is q^ < Q"^). Experimentally, the most important 
situation is when one of the photons is almost real q^ ~ 0. In this case, necessary nonperturbative 
information is accumulated in the same integral / (see eq.(|l])) that appears in the one-gluon- 
exchange diagram for the pion electromagnetic form factor |^, The value of I depends on 
the shape of the pion distribution amplitude (p-n-^x). In particular, using the asymptotic form 
(/?"*(a;) = QfnXX [|, ^ gives F"^.^o(Q^) = 47r/^/Q^ for the asymptotic behaviour [jT|. If one takes 
the Chernyak-Zhitnitsky form Lp^^{x) = 30/,ra;(l — a;)(l — 2x)'^ |^, the integral / increases by a 
sizable factor of 5/3, and this difference can be used for experimental discrimination between the 
two forms. 

An important point is that, unlike the case of the pion EM form factor, the pQCD hard 
scattering term for 77* — > vr^ (7 denoting a real photon) has zeroth order in the QCD coupling 
constant a^, i.e., the asymptotically leading term has no suppression. The situation is similar 
to that in deep inelastic scattering. Hence, we have good reasons to expect that pQCD for 
F^^*Tj.o{Q^) may work at accessible Q^. Of course, the asymptotic l/Q^-behaviour cannot be true 
in the low-Q^ region, since the = limit of F^^*^o(Q^) is known to be finite and normalized 
by the 77 decay rate. Using PCAC and ABJ anomaly 0, one can calculate F^^*t^o{0) 

theoretically: F^^.^o(O) = I/tt/tt. It is natural to expect that a complete QCD result does not 
strongly deviate from a simple interpolation vr/^F^^.^o((5^) = 1/(1 + (5^/4vr^/^) between the 

= value and the large-Q^ asymptotics^. This interpolation implies the asymptotic form of 
the distribution amplitude for the large-Q^ limit and agrees with CELLO experimental data [^]. 
It was also claimed []T0| that the new CLEO data available up to 8GeV^ also agree with the 



interpolation formula. This provides a strong evidence that the pion distribution amplitude is 
rather close to its asymptotic form. Because of the far-reaching consequences of this conclusion, 
it is desirable to have a direct QCD calculation of the 77* 7r° form factor in the intermediate 
region of moderately large momentum transfers Q^>lGeV^. Such an approach is provided by 
QCD sum rules. As we will see below, the QCD sum rules also allow one to calculate F^^*^o(Q^) 
for large without any assumptions about the shape of the pion distribution amplitude. In fact. 



^In particular, such an interpolation agrees with the results of a constituent quark model calculation jsj 
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the QCD sum rule for F^^*^o{Q^) can be used to get information about 



2. Definitions. 

The 7*7* — > 7r° transition form factor F^.^.^o(g^, Q"^) can be defined through the matrix element 
J {7T,P \T{J^{X) J,(0)} |0)e-^''i^d^X = aV2e^,^fsq^ql F^*^*^. [q\Q^) , (2) 

where a = /An is the fine structure constant, = e {^u'-^^u — \d'-)^d^ is the electromagnetic 

current of the light quarks and |vr, P) is a vr" state with the 4-momentum p. To incorporate QCD 
sum rules [|ll|, we consider a three-point correlation function 

^«^.(gi,g2) = ^/ (0|t{j^(F)J^(X) J,(0)} |0)e-^''^^e^^^rf^Xrf^F, (3) 

(cf. |]12[) containing the axial current = -^{u'^^'^aU — d^y^'jadj serving as a field with a non- 
zero projection onto the neutral pion state: (0|j^(0)|7r°, P) = —if^Pa- The three-point amplitude 
^afMuiqi, q2) has a pole for = ml: 

i.e., the Lorentz structure of the pion contribution is Pa^^MuapQiQi^ and the spectral density of the 
dispersion relation 

^ (p\ q\ Q^) =- r ^^^'"^''^'^ ds + ''subtractions^ (5) 

^ ^ TT JO S — p'^ 

for the relevant invariant amplitude can be written as 



P 



(^s,q'^,Q'^^ = TifT,5{s — m^)F^*^*7r° (q'^,Q'^^ + ''higher states". (6) 



The higher states include Ai and higher broad pseudovector resonances. Due to asymptotic 
freedom, their sum, for large s, rapidly approaches the pQCD spectral density p^^{s, g^, Q"^). The 
simplest model is to approximate all the higher states, including the Ai, by the perturbative 
contribution: 

p-"'^ {s, q\ Q') = 7rfJ{s)F,,,,^o (q\ Q^) + e{s - s,)p^^(s, q\ Q') (7) 

where the parameter So is the effective threshold for higher states. To suppress the higher states 
by an exponential weight exp[— s/M^], we apply the SVZ-Borel transformation ||11|| : 



5( V - M^) g^ Q') ^ $(M^ g^ Q^) = e-^/^^V(^, g^ Q') ds, (8) 
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Figure 1: a) Leading-order pQCD contribution, b) Triangle diagram, c) General structure of 
bilocal contributions. 

which, moreover, produces a factorially improved OPE power series for large M^: l/(— p^)^ 
(l/M2)^/(Ar- 1)!. 

To construct a QCD sum rule, one should calculate the three-point function g^, Q^) and 

then its SVZ-transform $(M^,g^,Q^) as a power expansion in for large M^. However, a 

particular form of the original ^-expansion depends on the values of the photon virtualities 

and Q"^. 
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-p^ > fi'^ where /i is some scale 



3. QCD sum rules for large q^. 

The simplest case is when both virtualities are large: Q 
/i^ ~ 1 GeV^ above which one can rely on asymptotic freedom of QCD. Then all contributions 
which have a power behaviour (1/M^)^ correspond to a situation with all three currents close to 
each other: all the intervals X^, Y^, {X — F)^ are small. 

The starting point of the OPE is the perturbative triangle graph (Fig.p^). Using Feynman 
parameterization and performing simple integrations we get: 



[sxx + xQ"^ + xq^Y 



dx . 



(9) 



The variable x is the light-cone fraction of the total pion momentum p carried by one of the 
quarks. Adding the condensate corrections, we obtain the following QCD sum rule: 



64 3 . ,2 1 

243 VM4 



= 2 Tdse-'/^'" f' ^^(^Ql±^SD^dx 
Jo Jo [sxx + xQ"^ + xq'^Y 

+-(^GG) ( ^ I ^ -) 

9 W ^ I 2M2Q2 ^ 2M2g2 Q2^2 ) 
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Q2g4 Q4^2 



(10) 
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It is valid when both virtualities of the photons are large. In this region, the perturbative QCD 
approach is also expected to work. This expectation is completely supported by our sum rule. 
Indeed, neglecting the sxx-term compared to xQ^ + xq^ and keeping only the leading 0(1/(5^) 
and 0{l/q^) terms in the condensates, we can write eq.(|TUp as 



)xx 



-^,(^GGmx) + 5{x)] 



+ +^(^)] +2[5'(a:) + . (11) 

The expression in curly brackets coincides with the QCD sum rule for the pion distribution am- 
plitude fn^n{x) (see, e.g., [|l3l). Hence, when g^, the smaller of two photon virtualities is large, 
the QCD sum rule (p!OD exactly reproduces the pQCD result (0). 



4. Operator product expansion for small q^. 

One may be tempted to get a QCD sum rule for the integral / by taking = in eq.(0). Such 
an attempt, however, fails immediately because of the power singularities l/q"^, l/g^, etc. in the 
condensate terms. It is easy to see that these singularities are produced by the 6{x) and 6'{x) 
terms in eq.(ll). In fact, it is precisely these terms that generate the two- hump form for ip-j^lx) 
in the CZ-approach p. Higher condensates would produce even more singular S^'^\x) terms. As 
shown in ref.[jl3[, the 5^^\x) terms result from the Taylor expansion of nonlocal condensates like 
(g(O)g(Z)). Modelling nonlocal condensates by functions decreasing at large (— Z^), i.e., assuming 
a finite correlation length ~ 1/yU for vacuum fluctuations, one obtains smooth curves instead of the 
singular 6{x) and 6'{x) contributions, and the result for ip.j^{x) is close to the asymptotic form [[l^] . 
Effectively, the correlation length provides an IR cut-off in the end-point regions x ~ 0,x ~ 1. 
Similarly, we expect that, when is too small to provide an appropriate IR cut-off in the sum 
rule ([T0|), such a cut-off is again generated by nonperturbative effects, i.e., that eventually 1/g^ 
is substituted for small g^ by something like l/m'^p. Below, we show that this is exactly what 
happens in the QCD sum rule framework. 

To illustrate the nature of the modifications required in the small-g^ region, it is instructive 
to analyze first the perturbative term. The latter, though finite for g^ = 0, contains contributions 
which are non- analytic at this point: 



$^V,Q^M2) = -^/ e-^/^^^{ 1 + 



M2 



where y = s/{s + Q^) and dots stand for terms analytic and vanishing for g^ = 0. The logarithms 



here are a typical example of mass singularities (see, e.g., [|1^, |T5| and, for QCD sum rule applica- 



5 





a) 



b) 



c) 



Figure 2: Bilocal contributions with coefficient functions given hj a) single propagator; b) product 
of 3 propagators; c) product of 2 propagators. 



tions, refs. |16|, |18|): the singularities are due to the possibility of the long-distance propagation 
in the g-channel. In other words, when is small, there appears an additional possibility to get a 
power-behaved contribution from a configuration in which the large momentum flows from the Q- 
vertex into the p- vertex (or vice versa) without entering the g- vertex, and small momenta flowing 
through other parts of the diagrams induce singular contributions. In the coordinate represen- 
tation, such a configuration can be realized by keeping the electromagnetic current Jfi{X) of the 
low-virtuality photon far away from two other currents, which are still close to each other. The 
contribution generated in this regime can be extracted through an operator product expansion for 
the short-distance-separated currents: T{J{0)j^(Y)} ~ J2Ci(X)Oi. Diagrammatically, the situa- 
tion is similar to the pQCD limit g^, ^ — discussed above. The only difference is that we 
should consider now the limit —p^^Q'^ ^ q^. The result again can be written in a "parton" form 
(see Fig.lc): 



(13) 

where 4>^\{y} ■, q^) can be treated as distribution amplitudes of the gi-photon, with y's being the 
light-cone fractions of the momentum qi carried by the relevant partons (i.e., quark and gluonic 
G-fields present in O). The functions 4>^2\{y},q^) are related to the correlators ("bilocals", cf. 



m) 



n«(?i^ 



;*«i^(0|r{J^(X)C«(0)}|0)c/'X 



(14) 



of the J^(X)-current with composite operators 0^'^\Qi). Performing such a factorization for each 
diagram, we represent the amplitude as a sum of its purely short- distance {SD) and bilocal 
[B) parts. The S'-D-part (which is defined as the difference between the original unfactorized 
expression and the perturbative version of its S-part) is regular in the — * limit and can be 
treated perturbatively. On the other hand, the low-g^ behaviour of the ^-correlators n(gi) cannot 
be directly calculated in perturbation theory. 
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5. Structure of bilocal contributions. 



In the simplest case, the amphtude T^^\{yqi}; q2,p) in the bilocal term is given by a single quark 
propagator (Fig. 2a): 

T{yqi,yqi;q2,p) ^ 1 — = — — ^ 7^, (15) 

[p -yqi) yp + vvq - vQ 

accompanied by two-body distribution amplitudes (f)iy\y, q"^), with yqi and yqi being the momenta 
carried by the quarks (y = 1 — y). The |/"-moments of (t>^^\y ■, q"^) are given by bilocals n«(g2) 
involving composite operators O^*-* with n covariant derivatives. Note, that it is legitimate to keep 
the g^-term in eq.(|T^) when substituting it into eq. (P^) : all the {q^/Q"^)^ and {q^/p^)^ power 
corrections are exactly reproduced there due to a phenomenon analogous to the ^-scaling |20| in 
deep inelastic scattering. Applying the SVZ-Borel transformation, we obtain the result 

$f (g^ Q\ M^) ~ ^ <P,{y, q')ey'^' l'^\-y<^' ly^' dy, (16) 

which has the structure of eq. (P^ : one should just take yQ"^ /y = s. 

In perturbation theory, the amplitudes (l)^\y,q'^) have logarithmic non-analytic behaviour for 
= 0. Eq.(p!2|) indicates that, for the triangle graph, there are only two independent sources of 
logarithmic singularities. They correspond to two-body operators of leading and next-to-leading 
twist in the OPE for T{J{0)j^(Y)}. The logg^-terms reflect the fact that the lowest singularity in 
$^^(g^, Q^, M^) for the g-channel corresponds to threshold of the gg-pair production which, for 
massless quarks, starts at zero. However, this is true only in perturbation theory. For hadrons, 
the first singularity in the g-channel is located at the tttt threshold, with the p-resonance being 
the most prominent feature of the physical spectral density for the correlators 11^*) (g^) . in other 
words, the 1/g^ and 1/g'^ terms in the condensates and logarithmic terms in the perturbative 
contribution correspond to the operator product expansion for the correlators 11^*^ (g^), which is 
only valid in the large-g^ region. To get nW(g2) for small g^, one can use this large-g^ information 
to construct a model for the physical spectral density (Jn\t) and then calculate 11^*^ (g^) from the 
dispersion relation 

1 a^^\t)dt 
n Jo t + q^ 

To this end, we use the usual "first resonance plus continuum" model 

^^it) = g^^Kt - nil) + e{t > .p)cr«^^(t). (18) 

Practically, this means that we modify the original spectral density in the region t < Sp hj 
subtracting all the terms of the OPE for an{t) in this region and replace them by the p-meson 
contribution g^^5{t — m^). In particular, this subtraction eliminates 1/g^ and 1/g^ singularities 
corresponding to the condensate 5{t)- and 5'(t)-contributions into cr„(t). For the perturbative 



n<-.tf) = if™ (17) 
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contribution, the subtraction procedure removes the logg^ terms. If cr^^{t) ~ t (this produces the 
log contribution) , the "correction term" is given by 

tdt g^^^ 2, (sn + q^\ g^^^ , , 

r + g"^ + \ h ) Q' + 

As a result, the — logg^ term cancels the first non-analytic term in the three-point function, and 
effectively one gets log Sp instead of log in the small-g^ region. In the large-g^ region, where the 
original OPE must be valid, the correction terms should disappear. Requiring that they vanish 
there faster than the contribution which they are correcting (i.e., faster than 1/g^) we arrive at 
the relation g'^^^ = In a similar way, if a^'^{t) ~ (this gives the g^logg^ term) one gets 

the relation g^"^^ = —3^/3. 

Imposing a universal n-independent prescription is equivalent to the assumption that the y- 
dependence of the p-meson distribution amplitudes coincides with that of the perturbative corre- 
lators and, furthermore, that the p-meson contribution is dual to the quark one, with the standard 
duality interval 3p ~ 1.5 GeV^ obtained in ref. |TT|. In principle, one can use more elaborate models 
for the distribution amplitudes (p^p\y). As far as the requirement of smallness of the additional 
terms in the large-g^ region is fulfilled, our results do not show strong sensitivity to particular 
forms of (p^p\y). However, using the local duality approximation for ^^p\y) we are able to present 
our results in a more compact and explicit form. 

After the modifications described above, the contribution of the triangle diagram converts into 

$nfM2 0^ O'^) = ^ r d3e~''^^^ — ll + e^'^"/*^' 



M2 y M2 " 2(m2 + g2) 



+ r M2 ' 2 3(m2 + g2) 



(20) 



(here y = s/(s + Q^)). As promised, in this expression, we can take the limit g2 ^ without 
encountering any non-analytic terms. Note, that the modified versions of g2 log g2 g^^^d g^ log g2 
terms do not vanish in the g2 — > limit. Finally, using the formula 

= 1 e-'/^'\g'{0)6{3)+g{0)6'{3)+g"{3))— (21) 

we can rewrite the g2 — > limit of eq. (^OD in the canonical form (H) and determine the relevant 
spectral density po{3,Q^). 

Most of the singular condensate contributions of the original sum rule (0) can be interpreted 
in terms of the bilocals corresponding to the simplest, one-propagator coefficient function. As a 
result, they are subtracted by the procedure described above. To factorize the gluon condensate 
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contribution, we used a technique similar to that developed for the perturbative term. For the 
diagrams with the quark condensate, the factorization in most cases is trivial. However, there are 
also terms with the coefficient function formed by three propagators (see Fig.26). In this case, the 
long-distance contribution is described by the photon distribution amplitude (p^{y,q'^) related to 
the O ~ g75[7a, 7/3]D . . . Dq operators. The OPE for such a "non-diagonal" correlator (cf. pl[] ) 
starts with the term proportional to the quark condensate: 4>^{y,q'^) ~ ■^{qq)[^{y) + ^{^ — v)], 
strongly peaked at the end-points. However, incorporating the nonlocal condensates to model the 
higher-dimension contributions and employing a novel technique applied earlier to a similar 
non-diagonal correlator, we found that the distribution amplitude of the lowest-state (i.e., p- 
meson) is rather narrow. Neglecting the higher-state contributions (whose distribution amplitudes 
have oscillatory form) we obtain in the q^ = limit 

*"(^^ Q') = ^T^Af£ r ^ds r ^dy, (22) 

^ ' 27m2Q4]Vf6 Jo A/(s+Q2) y^ ^ ' 

where (p'^{y) is the normalized distribution amplitude (its zeroth y-moment equals 1), which we 
model by (p'^{y) = 6?/(l — y). Transforming $"^(M^,Q^) to the canonical form by using eq.(|21|) 
gives the spectral density p^{s,Q'^). Note that, for our model, p^{s,Q'^) contains the (l/s)+ 
distribution. 

The OPE for the three-propagator coefficient function also produces operators like q . . . 7^Z^'^g. 
Naively, one would expect that such operators vanish due to equations of motion. However, when 
inserted in a correlator, they produce the so-called contact terms [^]. In our case, the contact 
terms give 



27Q6M6 Jo 



ln^-2- 



sds. (23) 



Again, the spectral density p'-^{s,Q^), obtained after applying eq. (|21|) , contains the (l/s)+ distri- 
bution. 

Finally, there are also configurations with the coefficient functions given by two propagators 
(see Fig. 2c) which correspond to three-body qGq-tjpe distribution amplitudes. Their contribution 
was found to be small and, to simplify the presentation, we will not include them here. 



6. Sum rule for F^*^*^o(g^, Q^) in the q'^ = limit. 

Since all the contributions, which were singular or non-analytic in the small-g^ limit of the original 
sum rule, are now properly modified, we can take the limit q^ = and write down our QCD sum 
rule for the 77* —>■ 7r° form factor: 
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+ 2 
+ 



g2 _6s + 3sVQ^ fs 



{s + Q 



\2\4 




-s/M 



-s/M 



64 3 ,2 f 1 12 



+ 



■so 



-s/M2 



/o 
4 



^5^ + 3gg^ + 4Q^ 



P L 
1 



2 Q^ds 








(is 








- 2 






(is 





log— -3 + 









s + 

s2 + 3sQ2 + 6Q4 



(s + Q 



2\3 



s + A^ 



(is 



(24) 



The sum rule must be taken in the limit A^ — of the parameter specifying the regularization 
which we used to calculate the integrals with the (l/s)+ distribution. Furthermore, this sum rule 
implies that the continuum is modeled by an effective spectral density p^-^-^(s, Q^) rather than 
by p^"^(s,Q^), with p'^ff{s,Q'^) including all the spectral densities which are nonzero for s > 0, 
i.e., po('5, Q^), P^is, Q^), p'"(s, Q^) and also an analogous contribution from the gluon condensate 
term. 

1.5 GeV^ 
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Using the standard values for the condensates and p-meson duality interval Sp 
we studied the stability of the sum rule with respect to variations of the SVZ-Borel parameter 
in the region > 0.6 GeV^. Good stability was observed not only for the "canonical" 
value Sq = 47r^/^ ^ 0.7 GeV"^ but also for smaller values of Sq, even as small as O.AGeV"^. Since 
our results are sensitive to the so-value, we incorporated a more detailed model for the spectral 
density, treating the Ai-meson as a separate resonance at s = 1.1 GeV"^, with the continuum 
starting at some larger value sa- The results obtained in this way have good M^-stability and, 
for < 1.2GeV'^, show no significant dependence on sa- Numerically, they practically coincide 
with the results obtained from the sum rule (^) for sq = 0.1 GeV^. 

In Fig.^], we present a curve for (5^-F^^.7ro(Q^)/47r/^ calculated from eq.(p4D for sq = 0.1 GeV^ 
and = O.'&GeV'^. It is rather close to the curve corresponding to the Brodsky-Lepage interpo- 
lation formula nf-KF^^tT^o^Q^) = 1/(1 + Q'^ /i-n'^fl) and to that based on the p-pole approximation 
7r/^F(Q^) = 1/(1 + Q"^ /rn?p). It should be noted, however, that the closeness of our results to the 
p-pole behaviour in the Q^-channel has nothing to do with the explicit use of the p-contributions 
in our models for the correlators in the g^-channel: the Q^-dependence of the p-pole type emerges 
due to the fact that the pion duality interval sq ~ 0.1 GeV^ is numerically close to 0.6 GeV^. 

For Q"^ < ?)GeV'^, our curve goes slightly above those based on the p-pole dominance and 
BL- interpolation (which are close to the data |^). This overshooting is a consequence of our 
assumption that can be treated as a large variable: in some terms, l/Q^ serves as an expan- 
sion parameter. Such an approximation for these terms is invalid for small and appreciably 
overestimates them for Q"^ ~ IGeV"^ producing enlarged values for F^^.^o(Q^). 

In the region Q"^ > 3GeV^, our curve for Q^F^^*^o(Q^) is practically constant, supporting the 
pQCD expectation (|T]). The absolute magnitude of our prediction gives / ~ 2.4 for the /-integral. 
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4nf 




Q' (GeV') 

Figure 3: Combination Q^Fj^*^o{Q^)/4:iTf.„ as calculated from the QCD sum rule (solid line), 
p-pole model (dashed line) and Brodsky-Lepage interpolation (dash-dotted line). 



Of course, this value has some uncertainty: it will drift if we change our models for the photon 
distribution amplitudes (bilocals). The strongest sensitivity is to the choice of <^T(|/) in the tensor 



contribution (p2D . However, even rather drastic changes in the form of (Pp{y) do not increase our 
result for / by more than 20%. The basic reason for this stability is that the potentially large 
factor from the relevant contribution in the original sum rule ([To|) is substituted in (^2]) by 
a rather small (and non-adjustible) factor l/m'^p- 

Comparing the value I = 2.4 with J*^* = 3 and J*"^ = 5, we conclude that our result favours a 
pion distribution amplitude which is narrower than the asymptotic form. Parametrizing the width 
of (pn{x) by a simple model (pn{x) ~ [x{l — x)]"', we get that / = 2.4 corresponds to n = 2.5. The 
second moment (^^) (^ is the relative fraction ^ = x — x) for such a function is 0.125. This low 
value (recall that {^Y' = 0.2 while (^^^cz ^ g ^g^j 

agrees, however, with the lattice calculation 
23 and also with the recent result obtained from the analysis of a non-diagonal correlator. 



7. Conclusions. 

Thus, the QCD sum rules support the expectation that the Q^-dependence of the transition form 
factor F^^*T^a{Q'^) is rather close to a simple interpolation between the = value (fixed by the 
ABJ anomaly) and the large-Q^ pQCD behaviour F{Q'^) ~ Q~^. Moreover, the QCD sum rule 
approach enables us to calculate the absolute normalization of the term. The value produced 
by the QCD sum rule is close to that corresponding to the asymptotic form V5"'*(x) = QfT,x{l — x) 
of the pion distribution amplitude. Our curve for F^^tT^oiQ"^) is also in satisfactory agreement with 
the CELLO data ||^ and in good agreement with preliminary high-Q^ results from CLEO [pil|] . 
Hence, there is a very solid evidence, both theoretical and experimental, that '{>-,t{x) is a rather 
narrow function. 
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